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ABSTRACT 


Unlike an earlier shock spectra approach, the generalization 
presented in this paper permits an accurate elastic interaction between 
the spacecraft and launch vehicle to obtain accurate bounds on the 
spacecraft response and structural loads. In addition, the modal 
response from a previous launch vehicle transient analysis — with or 
without a dummy spacecraft — is exploited to define a modal impulse as 
a simple idealization of the actual forcing function. The idealized 
modal forcing function is then used to derive explicit expressions for 
an estimate of the bound on the spacecraft structural response and 
forces . 

Greater accuracy is achieved with the present method over the 
earlier shock spectra, while saving much computational effort over the 
transient analysis. 
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SECTION I 


OVERVIEW 

A. BACKGROUND 

1. Transient Analysis 

The determination of the internal loads in a spacecraft structure 
subjected to a dynamic launch environment is an iterative process aimed 
at the ultimate sizing of the structure for its design. The process 
is iterative because the internal loads in the spacecraft are dependent 
upon the spacecraft mass and stiffness distribution, as well as the 
mass and stiffness distribution interaction with the launch vehicle. 

Some of the key parameters that dominate the loads analysis process 
are timeliness of the sizing, cost, weight limitations, and low 
sensitivity to iterative changes in the design. Because the loads are 
time-dependent functions, a complete transient analysis aimed at deter- 
mining responses to given time-dependent forcing functions is widely 
used throughout the industry to determine loads in the launch vehicle 
and in the spacecraft. The transient analysis is a simulation method 
by which the loads time histories must be computed in their entirety 
before peak values can be identified for establishing design loads. 

Although a transient analysis tends to result in lightweight 
designs when properly used, it is usually costly and time consuming, 
especially when several organizations are involved. No better method 
is presently known when dealing with a launch vehicle, and it is not 
suggested here that this process be changed for the launch vehicles. An 
alternate method is presented, however, for the analysis of the 
spacecraft . 

2. Alternates to Transient Analysis 

In the analysis of an aerospace structure, at least two organi- 
zations are involved: one organization responsible for the design and 
analysis of the launch vehicle, and another for the spacecraft. For 
the spacecraft-developing organization, it is desirable to perform 
frequent but inexpensive loads analyses to account for frequent design 
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changes in the spacecraft. Ideally, these analyses would include only 
the spacecraft, without requiring launch vehicle modeling and without 
requiring activity at the launch vehicle agency. For this reason, 
several alternatives have been suggested, and References 1-1 and 1-2 
are recent examples. Like the method presented here, the methods given 
in References 1-1 and 1-2 apply only when the loads in the spacecraft 
resuit from time-dependent forcing functions applied on the launch 
vehicJe, and bypass most of the costly time consuming steps of a com- 
plete transient solution. 

3. Reduced Transient Analysis Method 

The method of Reference 1-1 is a reduced transient analysis. 
Accelerations at the spacecraft-launch vehicle interface are obtained 
by a single transient analysis of the launch vehicle with a rigid mass. 
From this single launch vehicle analysis, interface accelerations 
can then be corrected to account for the elasticity of the spacecraft 
and used as input to the spacecraft. The correction procedure requires 
knowledge of the characteristics of the elastic spacecraft and of the 
rigid mass/elastic launch vehicle, but does not require the actual 
integration of the spacecraft elastic model with the launch vehicle 
elastic model, resulting in a significant cost reduction. From there 
on, the procedure follows the usual transient analysis. 

4 . Shock Spectra Method 

Like the approach presented here, the method of Reference 1-2 is 
a generalization of the traditional shock spectra concept. In this 
concept, shown in Figure 1-1, a shock spectrum S(w, £) of a function 
Xg(t) is the largest peak response of a single-degree-of- freedom 
oscillator subjected to the base forcing function, Xg(t). The implicit 
assumption is made that the input forcing function, Xp(t), is not 
affected by the presence of the escalator. This is true only if the 
mass of the oscillator is infinitesimal. 

The most significant feature of Reference 1-2 was the introduction 
of the relative impedance between the spacecraft and the launch vehicle, 
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Figure 1-1. Shock Spectrum S(u>, 4) of a Function x ( ,( t ) 
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i.e., between the oscillator and the base, in the form of a reduction 
factor, C, which accounts for the effect of the oscillator on the inter- 
face acceleration. In addition, an envelope is constructed for the 
acceleration spectrum at the interface between the spacecraft and launch 
vehicle. Figure 1-2, which is used to determine the level of an equiva- 
lent sinusoidal forcing function to be applied at the interface for the 
computation of the spacecraft member loads. In this manner, spacecraft 
loads are computed with no integrated launch ide/spacecraft transient 
analysis. This method was applied to several spacecraft structures as 
reported in Reference 1-3, and was found to be much less costly to 
use, but produced somewhat heavier structures than a transient analysis. 

In the application of the early shock spectra method outlined 
above , a number of assumptions were introduced that made its successful 
usage highly dependent upon the ane^yst's intuitive ability to recognize 
launch vehicle modes that are important to retain, and those of insig- 
nificant effects that can be discarded. The method sensitivity to the 
analyst's judgment is highlighted in the tacit assumption that the 
spacecraft /launch vehicle interface response is dominated almost totally 
by one launch vehicle mode at any one frequency, and in the process by 
which the envelope of the shock spectra is constructed. 

5. The Generalized Shock Spectra Method 

The present approach is an extension of Reference 1-2, with the 
objective of eliminating the analyst's judgment, formalizing the 
derivation and computational steps of the approach, explicitly stating 
the necessary assumptions, relaxing unessential ones, and introducing 
several significant improvements. These will be discussed in subsequent 
sections . 

B. OUTLINE OF THE GENERALIZED SHOCK SPECTRA 

The following general observations are fundamental to understand- 
ing the rationale of the shock spectra approaches: 

(l) The general objective is to avoid the cost of a launch 
vehicle/new spacecraft overall transient analysis. 
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(2) There has been a previous dynamic analysis done for the 
launch vehicle with another or a dummy spacecraft, the 
results of which are used as inputs to the generalized 
shock spectra method. Note that the present method is 
derived for use of the modes of the launch vehicle loaded, 
by a rigid mass only, or not loaded. If such an analysis 
is not available, the proper data can still be recovered 
from the analysis of the launch vehicle loaded with an 
elastic spacecraft as shown in A;pendix C. 

(3) The structural analyst needs only to determine the worst 
case load maxima or bounds, rather than time histories, 
of the structural response. Since maxima or bounds are 
the objective of the shock spectra concept, a shock spectra 
approach is readily applicable to the structural design 
process. 

(4) Any loads analysis — transient, shock spectra, or other — 
incorporates two basic items: a model idealizing the 
dynamic environment, i.e., the forcing function, and another 
idealizing the composite structural system. The form of 
each model is influenced by the selected approach, while 

its complexity is constrained by cost and time. 

(5) The original launch vehicle forcing functions, usually 

unknown to the spacecraft analyst, may be idealized by a 

simpler form to obtain an explicit closed form solution. 

A complete definition of this idealization is possible from 
the previous analysis of item (l) above, and is done at 
the modal level . 

In view of these remarks, the generalized shock spectra method 
is highlighted as follows: 

First , the modal forcing function F^(t) corresponding to the 
modes of the launch vehicle loaded by a rigid mass at the spacecraft 

interface, and representing the modal contribution of an actu . flight 
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event, is modeled, regardless of its physical point of application, by 
an equivalent launch vehicle modal forcing function, F (t). Unlike 
the actual complex transiert force, the equivalent forcing function 
chosen is a simpler form of variation with time. Here, an impulse 
delta function F e ^(t) * with a yet-unknovn magnitude F , or 

equivalently a velocity with magnitude v Q£ , is chosen for convenience. 

The choice of a simplified forcing function as an impulse 
emphasizes the view that the shape of the response time history is of 
no consequence, and that only the peak or bound of the response is of 
interest. Therefore, any forcing function that would reproduce a 
response with the same maximum peak or bound as the actual forcing 
function is acceptable. Tne equivalency between the actual forcing 
function and the idealized one is established, not on the basis of 
producing identical response time histories, but on the basis of pro- 
ducing an identical peak of the shock spectra of each of the launch 
vehicle modal response q^t). Figure 1-3, derived from the previously 
performed launch vehicle/dummy spacecraft analysis. An alternate to the 
direct use of q^(t) is discussed later. 

The use of modal shock spectra, rather than the interface uegrees- 

of-freedom shock spectra is significant because it automatically 

accounts for the matching of all interface physical degrees of freedom, 

and allows one to determine the modal magnitude of the impulse F^ or 

velocity v . It is noteworthy that the above process of establishing 
u x» 

the equivalent, idealized forcing function requires knowledge of the 
modal properties of the launch vehicle with or without a payload. Such 
information is usually available from the launch vehicle organization. 
Also note that the rigid mass at the launch vehicle /spacecraft interface 
does not have to be the total mass of the spacecraft to be analyzed, 
but can have any value convenient for the transient analysis performed 
earlier on the launch vehicle. The present approach does correct for 
whatever mass value was used. 

Second , in considering the composite structural system of 
Figure 1-2, which consists of a spacecraft modeled by S-normal modes 
and a launch vehicle modeled by L-normal modes, there will be (S +' L) 
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modally coupled equations of motion. Unlike the transient analysis 
where the solution is expressed in the complete (S + L) space of modal 
coordinates and time, a bound on the complete solution is established 
herein by: 

(1) Idealizing the totality of (S + L) mathematical space 

of modal coordinates by an array of nested (S x L) mathe- 
matical subspaces, in each of which only one spacecraft 
mode is individually coupled with one launch vehicle mode. 

In this fashion, each spacecraft mode is coupled with 
L-launch vehicle modes one at a time. To derive a bound 
on the total solution in the original (S + L) mathematical 
space, an explicit solution in the form of spacecraft 
modal response time history is first derived for the pair 
of modes in a typical subspace. The explicit form of this 
solution is derived here, and is based on the idealized 
modal forcing function just discussed. Furthermore, 
because the spacecraft member loads are the objective of 
the analysis, and since these are proportional to the 
generalized modal displacements, the response quantities 
used here are generalized modal displacements, q^ (t ) , from 
which an expression of the bound of q^Ct) is then 
derived. 

( 2 ) Numerical computations are made first to establish the 
bounds, Q , on each of the (S x L) discrete modal responses. 

S X/ 

Each bound Q corresponds to one of the (S x L) subspaces. 

To account for unknown design tolerances and variations in 
the structural model idealization, worst cases can be 
provided by allowing realistic possible tuning between each 
spacecraft mode and its nearest launch vehicle mode, and by 
scaling the entire frequency spectrum of the launch vehicle 
with respect to that of the spacecraft. 

(3) Next, a bound on the total spacecraft modal response, Q s , 
is constructed by summation over all the discrete L-bounds 
for that spacecraft mode. The summation can be over 
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absolute values, or in a root-sum-square sense that can 
also be weighted to account for phasing. 

(U) Finally, spacecraft member loads are obtained by adding 
the contributions of all spacecraft modes, either in 
absolute value or in a root-sum-square sense. 

In the procedure outlined above, steps (l) and (2) derive expres- 
sions for bounds on the discrete response, Q^, in each of the (S x L) 
modal subspaces, while steps / 3) and (U) construct a bound on the 
complete solution in the total (S + L) space from the discrete 
bounds. In this manner, much computational effort is saved over the 
usual transient analysis. 

The significant improvements or differences from Reference 1-2 are: 

(1) Idealization of the actual forcing function by an impulse 
permits a simple evaluation of the shock spectra peaks. 

Other forms of forcing function idealizations could also 
be dealt with in the least square sense. 

(2) The shock spectra matching is applied to the modal — rather 
than the physical — degrees of freedom. This gives 
greater flexibility and more generality in the application 
of the procedure. 

(3) The modal displacement spectra, Q g ^» is derived for an 
arbitrary pair of coupled modes, a launch vehicle mode 

and a spacecraft mode, regardless of the proximity of their 
associated natural frequencies. Thus, the a priori selection 
of a specific launch vehicle mode on the basis of its rela- 
tive significance in pairing with a given spacecraft mode 
is not required. 

(It) The entire treatment is believed to lend itself more 

readily to an automated spacecraft loads analysis that is 
substantially less dependent upon the analyst’s intuitive 
ability. 


ORIGIN ALP^EE 

op POOR QXJALm 


C. BASIC ASSUMPTIONS 

Other than the usual assumption of linearity of structural behavior 
for both the spacecraft and launch vehicle, the following assumptions 
ire made here: 

(1) The mass (and inertia) ratio between the spacecraft and 
the launch vehicle is small, but not necessarily 
infinitesimal. 

(2) The interface between the spacecraft and the launch vehicle 
is assumed statically determinate. 

(3) No external forces are acting directly upon the spacecraft 
other than through the common launch vehicle/spacecraft- 
interface . 

(b) A transient analysis of the launch vehicle with or without 
a dummy spacecraft is available. 

Assumptions 1 and 2 are not overly restrictive, and could be 
relaxed without much difficulty. Removing assumption 2 would require 
introducing more complex basic equations. Assumption 2 can also be dealt 
with by introducing a portion of the launch vehicle interface in the 
spacecraft ~*.del. Other lesser assumptions are also made, and will be 
stated in subsequent sections as they arise. 
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SECTION II 
BASIC EQUATIONS 

A derivation of the equations of motion of a damped linear 
structural system under the action of prescribed forces {F (t ) } and 

H *3 

prescribed rigid body accelerations {R(t}} can be found in the litera- 
ture (Reference 2-1). For each normal mode n = 1, 2, •••, N, the 
equations may be written in the following form: 


M (q + 2£ u q + w q ) = <4> >{F.} - <M >{R} 

nn ^n n n^n n hi nj j nR 


( 1 ) 


where , 


R = prescribed rigid body base motion 

= generalized modal coordinates that are related to the 

physical degrees-of-freedom, (u), by the normal modes, , 
such that {u^} = 


{ 4> } = normal mode shape associated with q ; it is designated by 
n "th ^ 

(<J> } when evaluated at the j physical degree-of- freedom 

J ^ 

£ = percent of critical modal damping for the n^* 1 mode 


a)^ = natural frequency of the system, associated with q^, {0^} , 

and 5 


<M _> = mass coupling between elastic and rigid body modes 
nn 


M = elements of the diagonal modal mass matrix whose normalization 
nn 

is subsequently chosen so that M * 1 for all modes 

nn 

{ F } ■ time variable vector of forces applied at an arbitrary location 
J 

j due to an arbitrary event 
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Equation (l) can be applied to a spacecraft mounted on a launch 
vehicle through a statically determinate interface as shown in Figure 1-2. 
Using subscripts s, i., and i to designate quantities associated, respec- 
tively, with the spacecraft, the launch vehicle, and their common 
interface i. Equation (l) can be separately written for the spacecraft 
and the launch vehicle by taking proper account of the reaction forces, 
{F^}, and accelerations, {R^}, at the interface. The results are: 

for the normal mode of the free-free launch vehicle 


H + + Vt = < ^j >{F J } ' < *t i >{F i } 


( 2 ) 


'or the s normal mode of the spacecraft, cantilevered from i. 


q + 2C u q +uq = -<m . >{R.} 
s s s s s s si 1 


(3) 


In Equations (2) and (3), the interface acceleration, {R^}, is 
obtained by superposition of the contributions of all launch vehicle 
modes, rigid and elastic: 

(i i> - t+J'V * [♦»)<¥*> (**» 

The launch vehicle rigid body acceleration {q } above is obtained as a 

n 

special case of Equation (2) when ui^ = 0: 

= “ ^Ri^ F i } ^ 

Furthermore, the reaction forces, { F^} , on the statically determinate 
spacecraft-to-launch vehicle interface are obtained by summing the 
reactions for all spacecraft modes T: 

[m is Hq s } s ■ 1, 2, •••, T (6) 
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in which [m^J is the contribution of each mode s to the rigid mass 
matrix relative to the interface i. The following relationship holds 
for each mode s: 


[m ii ] * {m is }<m si > 


Rather than including all spacecraft modes, only a truncated or selected 
set of modes S < T are usually retained. Such truncation or selection 
may be done by including a residual rigid mass [M. defined by 

11 nr*D 

(Reference 2-2 ) : 


^ii^RES 



( 7 ) 


where [M^] is the total rigid body mass matrix when all T modes are 
considered; 


[M. .] = 

li 


Ks 1 |m si J 


Therefore, Equation (6) for (Fj, can in general, be replaced by: 





(m S .) 

11 




( 8 ) 


In Equation (8), the neglected modes nr^ represented by their 
equivalent rigid mass, regardless whether their frequencies are higher 
or lower than the range in question. In this manner one is able to 
consider the elastic mass contributions of as few or as many spacecraft 
modes as desirable, with the remaining modes contributing to the 
residual rigid mass. A greater number of elastic modes result in a 
smaller residual rigid mass matrix, which approaches zero when all 
elastic modes have been included. From a modeling point of view, such 
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a residual mass matrix can also have a vide range of applications in 
representing the mass characteristics of part or all of the spacecraft 
during the various cycles of loads analysis. In particular, the residual 
mass term allows the analyst to use directly the modes of the launch 
vehicle loaded by a rigid mass. 

For ^spacecraft and L-launch vehicle elastic modes. Equations (2) 
and (3), along with their companion expressions (4), (5), and (8) result 
in (S + L) coupled equations describing the motion of the composite 
spacecraft-launch vehicle system. Instead of the commonly used transient 
analysis for the solution of this system of equations, a generalization 
of the shock spectra concept will be developed next, with the aim of 
obtaining a bound on the solution rather than its detailed time history. 
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GENERALIZED MODAL DISPLACEMENT SPECTRA 


A. IDEALIZED EQUATIONS 

In a full transient analysis of the coupled (S + L) Equations (2) 
and (3), all modes are simultaneously considered, and the solution is 
sought as time histories of the response quantities of interest. Since 
bounds rather than time histories are needed for design, much computa- 
tional effort can be saved by introducing mathematical idealizations 
that allow a relatively simple estimate of the bounds. This is achieved 
here when the entire (P + L) mathematical space of modal coordinates is 
idealized by an array of nested (S x L) mathematical discrete subspaces. 
In each subspace, a spacecraft mode s is coupled with a launch vehicle 
mode £, thus giving rise to (S x L) possible pairings. A bound on the 
true solution in the original (S + L) mathematical space is established 
by summation over the individual discrete bounds in the (S x L) 
subspaces . 


Because the spacecraft member loads are the object of the analysis 
here, and since these are proportional to the generalized modal dis- 
placements q , bounds on the spacecraft member loads are expressed in 
s 

terms of the bounds, Q , of the generalized modal displacements for 
each subspace. .A typical generalized modal displacement bound, Q , 
is derived here for a typical subspace in which one spacecraft mode q^ 
is coupled with one launch vehicle mode q^. For this purpose. Equations 
(L), (5), and (8) are solved together to yield: 


<V ■ (w ♦ ihiiKiO'Hi'h.ji'V + '»»> 


<V - ( [rl * KiAiO'HKiitV'V + + '"is 1 0 


( 9 . 




where 


Ki* - *VrES + Kil 

^iJ * [♦ 1 h H*£rW h1 ) 

['t'ij] * ^iRH^ R ][<l> R j] 

Tf the ratio of the Euclidian norm^ of the rigid mass matrix, [M^] 
defined above for the spacecraft, to the norm of the total rigid launch 
vehicle mass matrix, [Mpp] » is small, a simplifying approximation can 
be introduced in Equation (9), in which 

do) 

■ KiHV^V 4 ° 

This results in the approximate set of equations for and { } : 

( V - [tijHFj) * H u > i t - [* 11 ](» lB > % 

(li) 

{ V = [M ii K V * <m i s ) % 


which, when substituted into Equations (2) and (3), yields the following 
simplified version of the equations of motion for an arbitrary pair of 
coupled spacecraft and launch vehicle modes: 


^As a measure of the magnitude of a matrix A, the Euclidian norm is 
defined in analogy with the length of a vector by 


N(A) 


where a 


ij 


are the elements of A. 
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( 12 ) 


Note that removing the assumption that lead to Equation (10) will result 
only in redefining the terms of Equations (11 ) and (12). Also, although 
similar to Equation (C-5) of Appendix C, Equation (12 ) holds for only 
one launch vehicle mode and one spacecraft mode at a time, unlike 
Equation (C-5), which holds for all modes. 


In Equation (12), 


<*U >lM l l 1, *U ) * + < *li >l "i s ,<m si >U ill 1 


■ M sl * "st 


where 










M 


s2. 




”s£ * < ».i > [ M 0 ] ( » II! ) 






Then, Equation (12) can be rewritten in the following more convenient 
form: 


- 
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- 

(1 * "st 1 
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X 1 

» i. 
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i > 

X 2 


vvj 
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F,(t) 


(13) 
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in which the following additional definitions have been used: 


«. i 


m 


st 


si 1 + M 


s* 


1 / 1 + M 


si 


£ B — 

1 A + M 


si 


0) B U) 

2 s 


£ B F 
**2 


> (lU) 


F t (t) = 

4 * ' (r^vr) 



B. MODAL DISPLACEMENT SPECTRA 

The pair of modally coupled equations in Equation (13) describe 
the motion in any of the idealized (S x L) subspaces where a spacecraft 
mode is coupled with a launch vehicle mode. A bound on the spacecraft 
generalized modal displacement in such subspace is derived in this 
section from an explicit time-dependent solution of Equation (13), in 
which the actual launch vehicle modal forcing function F^(t) is modeled 
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by an equivalent forcing function having a simpler form of variation 
with time. The selected idealization of the forcing function is an 
impulse delta function having a magnitude F^, or alternatively, an 
initial velocity with a magnitude v^^. 

The equivalency between the actual forcing function and the idealized 
one is established on the basis of producing identical bounds on tl e 
modal displacement. This requirement yields one aquation for the 
definition of the magnitude of the initial impulse or velocity. The 
choice of an initial impulse or velocity as an idealization of the 
actual forcing function is made here only for simplicity of the resulting 
solution. 1> fact, any other choice of time variation for the idealized 
forcin • function is acceptable, provided that the abo\c mentioned 
equivalency condition is enforced. For example, a sinusoidal idealiza- 
tion (as in Reference 1-2) in terms of two unknown parameters, a magni- 
tude and a frequency, could be used. By enforcing more than two condi- 
tion equations, one could dete"inine the two unknowns, magnitude and 
frequency, in the least square sense. Functions of higher number of 
known parameters could be dealt with in a similar fashion. 

An analog computer solution of Equation ( 13 ) is given in Appen- 
dix A for a modal delta function impulse of magnitude F^, associated 
with a launch vehicle mode i.. Also, details of the analytical solution 
to the same equations are given in Appendix B for the equivalent modal 
initial velocity v Qj; . The analytical expressions for x^(t) and Xg(t) 
are written as: 


3-5 







oWSWW'tf^J* 

or took aJ“* 


The objective is to determine the maximum value of the displacement 
Xg(t) relative to the "base motion" x^(t). The maximum value of 
x 2 (t) is called the "generalized" displacement shock spectra since 
the "base acceleration" x^Ct) is alloved to be affected by Xg(t), unlike 
the simple traditional shock spectra concept. To derive relatively 
simple expressions requiring minimal numerical calculation for the 
generalized modal displacement spectra, further mathematical simplifica- 
tions will be introduced in Equations (15) and (l 6). These simplifica- 
tions were tested and found accurate when the results belov were 
compared with the results of the analog solution of Appendix A. 

2 2 

First, since and £ n g are much smaller than unity, the 
approximation 


for R t 0: 


r = f = f = 
^nl n2 n 


5 n! + ^2 


.-^fa a + M 0 1 
w„R L o 2 J 


for R = 0: 


5 

s n 2 


is used for the damped natural frequencies u> ^ and u> g, so that: 


where 


“„1 4 “ 2 (P 1 - ° 2 > 


“n2 4 “l <P l * P 2> 


“2 P 1 ‘ “2 S 0 


A - i 2 


damped center frequency 


Afi / 2 

"i.iPo * Wo v 1 - 


<T2 2 2 


damped beat frequency 
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Note that the approximation of Equation (17) is not used for the 
exponential part of the solution. Also note An/2Q Q « pg/p^ 

Furthermore, according to Equation (l6), 8 « 0. In fact, 8 tends 

to be exceedingly small in the neighborhood of B * 1. Therefore, it 

8t 

is sufficient to retain only the first order approximation for e , 
i.e., e®* * 1 ♦ 8t, and correspondingly e - ®* ■ 1 - 8t. This approxi- 
mation, along with that of Equation (17) and some trigonometric trans- 
formations, lead to the following simplified expression for the relative 
displacement x^(t) for the spacecraft mode: 

x 2 (t) A - w i v 0 i,l D 1 ( T ) + D 2* t ^* t = “ 2 t ( 2 °) 


where 


X I 1 — 8t 

D (t) = — T — I 0T I e sin p t cos p t 

2 w 2 P 2 1 1 


D (t) = 

2wf 


-1 ^1 

_ 2 p “ 

2» 2 »i \ °2 


“ 8 x 

0x | e cos p^T sin pgT 


For R ^ 0: since p^ is smaller than p^, an accurate estimate of 

the maximum value of the relative displacement x^(t ) of Equation (20) 

is provided by the discrete maxima D of the simple sinusoid, where 

cos p^t and cos PgT are either zero or unity, (see Figure 3-1). When 

D (t) approaches its maximum m aY > D 2 (t) approaches zero, and vice 

versa. Therefore, the maximum relative displacement D is given by the 

max 

largest of D, or D_ : 

Imax 2max 


D 


lmax ^ 2 

2 "l P 2 


D 


2max _ 2 

2 “ 2°1 



-0T . 

e sm 


-St . 

e sxn p 2 t 


V ( 21 ) 
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where t is determined as the smallest root of the characteristic 
equations 


for D. 


lmax ‘ 


tan P 1 t - g 


P 1 ^ “ ^1^2* 0T ^ 

Q ■ ■ ■ 

1 + (p i /p 2> b (1 “ Bt) 


( 22 ) 


for D, 


2max' 


tan Pgt * -g- 


i 


[l - (pg/p^ ©t] 

1 ■* (pg/pj^) j (1 - Bx) 


] 


The definitions introduced in Equation ( lU ) and the equalities 
of Equations (21) and (22) are now combined to give the generalized 
modal displacement spectra, Q gfc : 



Vot^st D 

(l + M^) max 


(23) 


Equation (23) above, along with its companion expressions of Equa- 
tions (21) and (22) represent a generalization of the traditional 
shock spectra concept since these expressions account for interactions 
between the oscillator and its base, i.e., between the spacecraft and 
launch vehicle modes. They also represent an improvement over the 
results of Reference 1-2 in that they hold true for an arbitrary pair 
of coupled modes, a launch vehicle mode and a spacecraft mode, regard- 
less of the proximity of their respective natural frequencies . It is 
this fact that makes it possible in the present approach to include 
contributions of all launch vehicle and spacecraft modes in the com- 
putation of bounds on the modal displacements and subsequently on the 
member loads . It is noted here that Equation (21) suggests that the 
damping effect is through the "average" damping 0 of the spacecraft 
and the launch vehicle as defined in Equations (l6) rather than through 
each damping separately. This point is explained further in Appendix A. 
Although derived for R / 0, Equation (23) is still valid in the limit 
for R -*■ 0. This is shown in Appendix D, 
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The generalized shock spectra will be identical to the traditional 
shock spectra in the limit when, in Equation (23), the mass ratio, y fljt 
defined by Eq. (l4) for the oscillator equivalent spacecraft mode, is 
allowed to approach zero while the associated stiffness increases such 
that the natural frequency is held unaltered. This is the case of the 
shock spectrum of a rigidly-loaded launch vehicle mode. This remark will 
be useful subsequently, in evaluating the initial velocity v^^ of 
Equation (23). 

Generally, however, for each set of parameters R, u^, C 1 , 

***, etc, describing the characteristics of a spacecraft mode s, and a 
launch vehicle mode 1, Equations (21), (22), and (23) will yield an 
estimate of the upper bound of the corresponding generalized modal 
displacement spectrum, Q g ^, in terms of the as-yet unknown initial modal 
velocity, v^. Once v^ has been determined, the generalized modal 
displacement spectra, Q gjl , for s * 1, 2, ••• S, and l = 1, 2, ••• L 
can be combined to establish a bound for the member loads. 

C. DETERMINATION OF THE INITIAL MODAL VELOCITY v QJl 

An analysis of the launch vehicle with some representation of the 
spacecraft is usually done by the launch vehicle organization for the 
structural design of the launch vehicle, or part of it. To properly 
design the part of the launch vehicle at the interface with the space- 
craft, the reaction loads from the spacecraft to the launch vehicle must 
be introduced in the analysis . Since these loads are almost entirely 
due to the fundamental cantilevered modes of the spacecraft, a simple 
model of the spacecraft such as a two-mass model is adequate. Note 
that this simple spacecraft model is used for the design of the launch 
vehicle at its interface with the spacecraft, but is not intended to 
be used for the detail design of the spacecraft. However, this two- 
mass model is useful for checking the launch vehicle data and can also 
be used to determine the weighting factor defined later. 



The original launch vehicle forcing functions as used for the 

launch vehicle analysis are given by a set of components { (t ) } expressed 

in the physical coordinate system xyz and applied at given locations of 

the launch vehicle structure. The modal analysis done on the launch 

vehicle with any spacecraft or dummy provides, in effect, a new system 

of coordinates, i.e., the generalized coordinates q. c (t), that are used 

to express new components of the forcing functions as a set of generalized 

forces F (t). In effect, F (t) is another definition of the same forcing 
c c 

functions F.(t). Therefore, one can use thi3 set (F (t)} to represent 
J C 

the forcing functions without the need of knowing the physical components, 
nor the location of the original forcing functions. Note that this 
observation holds only because an analysis of the launch vehicle with 
a dummy spacecraft has already been done. 

Any farther change of modes, if necessary, will simply change the 
set f F q ( t ) } by a modal transformation. For example, the present method 
requires the use of the modes of the launch vehicle loaded only by a 
rigid mass. Therefore, the components (F £ (t)} of the forcing functions 
with respect to the rigid mass loaded launch vehicle modes are 

(F £ (t)} = [V T ]{F c (t)> (24) 

where [V] is the modal matrix that transforms the elastically loaded 
launch vehicle to the rigid-mass loaded one as shown in Appendix C. 

Similarly, the generalized coordinates (q^(t)} of the rigid mass 
loaded launch vehicle are obtained from the coordinates (q. c (t)} of the 
elastically loaded launch vehicle by Equation (C-18) or Equation (C-20) 
of Appendix C. 

The unknown initial modal velocity v can be determined from 

U A' 

the requirement that the "actual" forcing function and the "idealized" 
forcing function produce identical bounds for the modal displacement 
spectra. Information for the "actual" forcing function can be obtained 
either from an initial analysis done on a simplified model of the 
spacecraft-launch vehicle combination, or from recorded flight response 
data for another spacecraft. It is believed that the minimum required 
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infornation is three degrees of freedom at the spacecraft-launch vehicle 
interface. The knowledge of more degrees of freedom improves the 
definition of v^^. In case of an available initial analysis, the 
knowledge of the modal response-time histories for each mode leads to 
the best estimate of v . For the last case, the peak value of the 

OX. p 

traditional displacement shock spectra, D^, is calculated for the modal 
acceleration response q^ft) of the rigid mass loaded launch vehicle as 
shown by Equations (28) and (29)» 

p 

Finally, it should be noted that will need to be computed only 
once for each launch vehicle mode and each flight event of interest. 

This part of the computation will remain the same as long as the launch 
vehicle and the forcing functions are not changed. 

In the derivation done so far and in what follows, rigid mass 
loaded launch vehicle modes are used. These modes are obtained by the 
method of Appendix C, or directly from the launch vehicle analysis if 
the spacecraft was represented by a rigid mass, as may be the case 
for preliminary launch vehicle design. 

1. Estimate of v QK From the Modal Response 

Assuming now that the time histories q^(t), l » 1, 2, ••*, L for 
the rigidly-loaded launch vehicle are available from a previous launch 
vehicle analysis, and that (w^, £ ) representing the peak of the modal 
displacement spectrum D (u> 9 , £~) for each mode has been evaluated, then 
Vqs, can be determined by requiring that be equal to its counterpart, 
Q^, derived from Equation (23) for M . ■+ 0. This is illustrated in 

X. P S*r 

Figure 3-2. To derive Q £ for R # 0, the shock spectrum Q 9 is computed 
for each of the launch vehicle elastic modes l - 7, 8, L from Equa- 

tion (23) in the limit as |i + 0. Appendix D gives the proper values 

S 

of the modal shock spectra for the case when R = 0, that is, for 

£• = 1, • • • , 6. 


* 


m V 

ion 


Lim 


(D ) 

max 


(25) 
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SHOCK SPECTRA 



Figure 3-2. Modal Displacement Spectra Dg and Q g 
for a Launch Vehicle Mode l 



The quantity, Lim (D ), of Equation (25) is readily evaluated by 

ji -K) max 

si 

substituting the following parameters in Equations (21) and (22): as 
U a •* 0, and R 4 0 

8 X 


e - I (^r + c 2 ) 


6 " 2 (C 1 R ' V 


2n o 

= (R + 1)/(R - 1) 


*i - 1 (a * D /i - 5* 


p = — (R - l) 

2 2 n 


A - S 2 




X (26) 


As seen from Equation (26), the displacement shock spectrum of Equa- 
tion (25) is dependent upon the value of R = u> /w . The peak value of 
P P Id p 

Q^, labeled Q^, occurs for R -*■ R . Finding the exact value of R that 

maximizes is rather involved analytically. Only an approximate 
estimate will be given here. 

2 

Considering the special case when G and of Equation (26) are 
assumed negligibly small, we find that 


-[tan 1 (a /2g)/(a /"B ) ] 

Lim (D ) = K 

U +0 maX 2 f~2 ~2 

sP. u)„ a + 


(r = 1 , 2 ) (27) 


2 


r 


where 


I (R + 1) for r ■ 1 

r 1 (R - 1) for r ■ 2 

By inspection of Equation (27), it is concluded that the peak occurs 

for r * 2 and R 2 * 1. The approximation that 0 in Equation (26) is 

negligibly small, is of the same order as the approximation that the 

imaginary part of the solution of Equation (15) is negligible, both of 

which are acceptable. The exact solution was evaluated from the analog 

simulation of Appendix A, in which the maximum of D occurred within 

■p max 

better than 1 % from R * 1. For R * 1, Equations (25) and (26) yield 
the following expression for the initial modal velocity and for 4 0: 

Vo* = eu) i (c i + V D * ( V V {28) 


where e = 2.718 . 

p 

Instead of the modal displacement peak D s (io , £ ), the modal 
P ^ 

acceleration peak A fl ((o„, £„) may be used. At the peak of the displace- 

x, 2 2 2 P P 

ment, the two quantities are related by V 5 *^’ ~ A^(Ug, . Thus, 

the following alternate expression for v^ is used, noting that 


Vo* 


= e(5, + £) 


Vv 


o 


(29) 


for 4 0 and where £ is the damping used in calculating the accelera- 
tion shock spectra of q^(t). 



As shown in Appendix D, 

v ) = A P 

0*'w -*0 * 


for io *0 Equation (29) becomes 


2 . Estimate of v^ From Physical Interface Response 

In the absence of availability of q c (t), a set of initial modal 
velocities v can also be obtained in a manner similar to that of 
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Reference 3-1 • It Is assumed that t _ time histories from an analysis 
or flight data of an degrees of freedom at the launch vehicle/spacecraft 
interface are known for another spacecraft or a dummy one. The shock 
spectra (i * 1, 2, N^) for the Nq degrees of freedom are calculated, 
as typically represented by the solid line of Figure 3-3. 

Next it is assumed that the mode shapes <j> of the structure corre- 

n 

sponding to the available interface acceleration are also known. Then 
each mode is given an initial velocity v Qn and the total interface 
acceleration is calculated from: 


Mt) = 


-E 

n*l 


to v_ 4> e 
n On n 


-S id t 
n n 


(a sin u>,t + b cos ui t) (30) 
n d n d 


where 


= (1 - 2C? )/ A - S 2 


b 

n 


2 £ 

n 


h) 


d 


w 

n 



The value of v. can be determined by trial and error so that 
On 

the shock spectra of R^(t) above, represented by the dashed line in 
Figure 3-3, envelopes the shock spectra of the real response. Since 
only a small number of vehicle modes are significantly contributing to 
R^, one can choose a set of modes N < N to simplify the calculations. 
The important criteria is that the v and the mode retained give an 
envelope of the real shock spectra. 

Finally, the initial modal velocity v QJl for the rigid mass loaded 
launch vehicle is obtained by the transformation of Equation (31), 
where a bound for v^ will be obtained by taking root sum square 
values, since the sign of v ^ is not retained. 


i 
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( 31 ) 


Note that this method is only a fall-back method, and that the modal 
response method of the previous section is to be preferred. 

D. TUNING 

During the early stages of design, the spacecraft and launch 
vehicle modes and frequencies are obtained from analyses that usually 
contain a large degree of uncertainty. To account for such uncertainties, 
one may introduce an artificial tuning between the spacecraft and launch 
vehicle modes. Unlike other methods, the present approach makes artifi- 
cial tuning possible and easy to implement because Equation (23) is 
inexpensive to evaluate and is valid for any pair of spacecraft and launch 
vehicle modes, regardless of the proximity of their respective 
frequencies. 

Two forms of artificial tuning have been identified: global and 
local. In global tuning, the entire spectrum of launch vehicle fre- 
quencies is incrementally ocaled in either direction V’<* to the 
spacecraft frequency spectrum. For each increment, a global response 



is computed and used as a measure fcr determining the worst case for 
design purposes. In this scheme, tuning is achieved by finding the 
amount of relative scaling that maximizes Q. Clearly, limits on 
the allowable relative scaling must be selected in advance, and the 
search for the maximum Q conducted within these limits. 

In the local tuning, the response is maximized for each spacecraft 
mode, one at a time. This is achieved by allowing the nearest launch 
vehicle frequency to coincide with that of the spacecraft frequency 
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under consideration) (i.e., R -*■ l) provided that the two were 
originally separated by no more than a preselected amount. Other 
schemes for tuning can be also devised. 



SECTION IV 




C 

O* SPACECRAFT MEMBER LOADS, 


DISPLACEMENTS, AND ACCELERATION 


A. MEMBER LOADS 

When the time histories of the spacecraft generalized modal 

displacement q (t) are determined as a result of a transient analysis, 
s 

one can express the vector of the spacecraft modal member loads, 

(f ( t ) } , for each member A due to motion in mode s in the general 
as a 

form (Reference 2-1 ): 


{f 

as 


(t)}. = [C JU } 
A ab bs 


% {t) 


(32) 


where 


= 1 , 2 , 


S 


In addition, the total load vector due to contributions of all space- 
craft modes is 


S 

{r a (t) L= £<Ls ( ‘»A <33) 

S=1 

where 

^ab^A = ma ^ r ^ x of> *' orce coefficients, whose elements are the 
th 

a force component in member A due to a unit displace- 
ment in the b th degree of freedom 

(4^,,} = as before, spacecraft mode shape at the b^* 1 degree- 
of- freedom 

[ C , ] ( 4>, ) = is recognized as the vector of modal stresses for mode s. 

ab bs 

In the present approach, the time histories of the spacecraft 
generalized modal displacements q (t) in Equations (32) and (33) are 



not computed for the mtire (S + L) mat ' at i cal space. Instead, they 
are replaced by an estimate of their upper bound, Q , so that the 
corresponding bound on the spacecraft modal member loads, is 

written analogous to Equation (32) in the form: 

<F .A • Q s (3 


The bound, Q , for the spacecraft generalized modal displacement 
s 

is expressed in terms of the discrete modal displacement spectra, Q ^ 
of Equation (23) for each of the (S x L) discrete mathematical sub- 
spaces. Since each results from coupling between a spacecraft mode 
s and only one launch vehicle mode, and since a complete representation 
of the launch vehicle includes more than one mode, say L modes, con- 
tributions due to all L modes should be included. Clearly, because 
Equation (23) does not contain information regarding time-phasing among 
modes, only a bound can be computed, and is provided by summation over 
absolute values or in the root-sum-square sense. 

Thus, Q is estimated by 
s 

L 

«. ■ 52 K* 1 <35 

*=1 

or 



or 


Q_ 


52 [ vvv 


N 


*=i 


(37) 


f . 3 « 
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whore W £^“ a » w £^ is a lighting function that can be used to account for 
time phasing between the launch vehicle modes. This weighting function 
can be determined from the available transient response of the simplified 
two-mass model of the spacecraft or from matching of the shock spectra 
of the mass-loaded launch vehicle interface. 

Analogous to the calculation of Q , a bound on the magnitude of 

s 

the total vector of member loads (F } due to all spacecraft modes 

A A 

s * 1, 2, •••, S is found from: 


{F ). 
a A 


S 

■ £ 

S®1 

» 


(F } « 

a A 


{F } = 

a A 


[W (w ) F ]‘ 
s s as 


Member loads computed by Equations (35) and (38) are, of course, much more 
conservative than those computed by Equations (36), (37), (39), and 
(40). Such excessive conservatism may be unnecessary, especially if 
each Q g for a spacecraft mode s is artificially tuned to the nearest 
launch vehicle mode. 

B. DISPLACEMENTS AND ACCELERATIONS 

An estimate of the bound on the displacements and accelerations 
associated with any or all degrees-of-freedom used in the spacecraft 


model can be obtained in a manner similar to the member loads of 
Equation (3b). Once Q has been established according to Equations (35), 

3 

(36), or (37), a bound on the spacecraft modal relative displacements 
(Dfag), or absolute accelerations {A^ g }, are found, respectively from: 



} 




(U) 


<\a> 


'“s Q s 


( 1 . 2 ) 


Finally, a bound on the spacecraft total displacement and acceleration 
degrees of freedom, (D^) and (A^) , is written analogous to Equa- 
tions (38) through (bO). For the displacements. 


<V* 



(1*3) 


or 


<V 



(VO 


or 


'V 



H5) 


= < 


Similarly, for the accelerations. 


ORIGINAL PAGE $ 
OP POOR QUAlfl* 


or 


or 



(U6) 




(1*7) 
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S 

II 

A 


H >W > \s |2 


N 

S=1 

• 


( 1 * 8 ) 


C. FLOW DIAGRAM 

Figure U-l summarizes the different steps taken from the launch 
vehicle data to the computation of bounds for the member loads, dis- 
placements, and accelerations. 
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Figure U-l. Generalized Shock Spectra Method 








SECTION V 




i 


np I £^ AL >,AG ® ® 

OF POOR QUALITY 


CONCLUSIONS 


The primary drawback of the original shock spectra method 
(Reference 1-1) is that it can result in very conservative answers , 
and therefore overdesign. It is also very dependent upon the analyst's 
judgment. However, it has the advantage of low cost of analysis and 
yields bounds that can be readily used for design. The transient 
analysis, on the other hand, is accurate, but is very expensive and 
difficult to implement because of its sensitivity to design change. The 
generalized shock spectra approach described in this paper, although 
still an approximate procedure, combines the advantages of low cost 
while maintaining a very reasonable accuracy. In addition, unlike the 
method of Reference 1-1, it is substantially less dependent upon the 
analyst's intuitive judgment in pairing spacecraft and launch vehicle 
modes. Because of its low cost, the tuning effect of launch vehicle 
modes with spacecraft mode can be very easily explored and used for 
design to establish worst case. Finally, the procedure is currently 
being applied to the Galileo spacecraft loads analysis with very accept- 
able results. 
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APPENDIX A 

ANALOG SOLUTION OF THE EQUATIONS OF MOTION 
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The analog exact solution of Equation ( 13 ) is given here to verify 
the analytical solution of Equation (15) and Appendix B. The equations 
of motion for the free vibration cf the system in Equation (13), or 
that of Figure A-l, are 

(1 + VI ) + Mx^ + 2£ 1 w ] x 1 + <^ 1 x 1 * 0 (A-l) 

2 

Ux^ + ux 2 + 2 ^ 2 ^ 2*2 * V)U> 2 X 2 * ° ( A ~ 2 ) 


where 




ratio of a spacecraft mode effective mass m 2 to a launch vehicle 
mode corrected effective mass m^; also, p = p g ^ of Equation (13) 

corrected natural frequency of one launch vehicle mode 

cantilevered natural frequency of a spacecraft mode 

launch vehicle modal damping 

spacecraft modal damping 

mass 1 response (unit mass) 

mass 2 response (mass p) 


The solution of Equations (A-l) and (A-2) to an impulse force 
applied on mass 1 is sought. Equivalently, the initial conditions for 
Equations (A-l) and (A-2) are 


x.^0) = F q x 2 (0) = 0 

x^O) = 0 x 2 (0) = 0 


(A-3) 


where F^ is a constant. 

To apply the initial conditions of Equation (A-3), the right 
hand side of Equations (A-l) and (A-2) must be 


F x = (1 + p) F q S(t) 

F 2 = pFq 6 (t ) 


(A-U) 


A-2 


Figure A-2 shows the analog diagram corresponding to Equations (A-l) 
and (A-2) solved by applying initial conditions of Equation (A-3). Fig- 
ure A-3 shows the analog diagram corresponding to Equations (A-l) and 
(A-2) solved by applying the impulse forces of Equation (A-4). These 
two circuits have been shown to give identical response time histories. 

All the results shown here have been obtained using Figure A-3. 

The following parametric variations were made: 

(l) Frequency ratio: 



f was held at 50.00 Hz and f^ was varied from 40.00 Hz 
to 60.00 Hz. 

(2) Dampings: = 0.01, 0.02, 0.03 and = 0.01, 0.02, 

0.03 maintaining ^ = constant = 0.04. 

(3) Mass ratio: p = 0, 0.001, 0.01, 0.1, 1.0. 

Figure A- 4 shows an example of the time histories of the spring 
force per unit mass w^x^t) for mass 2 versus time for all the param- 
eters R, C , and p . *ne quantity of interest for structural loads is 
^ 2 2 

the largest peak value D of to.x^(t) of mass 2, which is also the 

2 max 2 2 

peak of the total acceleration a of mass 2. "’he peak value was read 
r m 

from figures similar to Figure A-4 . Figure A-5 is an example of plots 

2 

of the peak value D versus the frequency ratio R. The results for 

2 max 

various values of the mass ratio p and dampings £ and £ 0 , keeping 
* 0.04 are tabulated in Table A-l. 

The solid curve on Figure A-5 has boon obtained from a simplified 
Equation (23) where 0=0 and where has been dropped in the damped 
natural frequencies 


a 

m 


D 

max 


“l F 0 


h 


An + 4 b 


-a(23/Aft) 

e 


(A-5) 


A-3 



where 


A£2 ■ /p + (1 - R 2 ) 

2B «= ? X R + (1 ♦ y) C 2 


a ■ tan ^ (AW/28) 


R * oJ^/Wg 


In addition, 


. . . An 

sin p^t = sin — t 


AH 

tan — t 


. 2 AQ 

V 1 + tan — t 


has beer used to show that the results depend on the square root of the 
mass ratio p since An contains p. The conclusions are as follows: 

(1) For a given value of mass ratio p, the maximum acceleration 

occurs for almost equal frequencies as is evident 

from Figure A-5. 

(2) Within an error of 1 % or less, the peak acceleration 
depends on the sum £ + and not on each separate 
damping. 

(3) The solid curve on Figure A-5 shows that Equation (A-5) 
is an excelle c estimate of the peak value of the 
acceleratior . 


The reduction factor C as defined in Reference 1-2 can be obtained 
as follows: 


a (p # 0) 

c = — 

a rn (p = 0) 


(A— 6 ) 


from the analog data of Table A-l. 



It can also be calculated by (jJlbBVlTuting Equation ( A -')) into 
Equation (A— 6 ) . 


C a 


* V 

An 2 + l+B 2 


1-a* ( 2g/Afi ) 
e 


( A— 7 ) 


where An, 6, and a are defined as for Equation ( A— 5 ) • Equation ( A— 7 ] 
provides the frequency dependence for C. 

Figure A-6 shows a comparison between the analog data and Equa- 
tion (A-7) as a function of p for various dampings .nd keeping 
+ S 2 = 0.0U, and for R * 1. This figure shows that Equation (A-7) 
is e very good estimate of C and that C ir virtually independent of the 
makeup of the sum Only one curve is shown for C as calculated 

by Equation (A-7). This is true because the value of C, when and 
are varied while keeping = 0.0*4, gives a difference too small 

to be plotted. 




Figure A-l. Launch Vehicle - Spacecraft Modal Model 
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Figure A-2. Analog Schematic for Initial Velocity Solution 
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Figure A-3* Analog Schematic for Impulse Solution 
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Figure A-5. Generalized Shock Spectra 
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Figure A-6. Reduction Coefficient C 
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Analytical solution is derived here in detail for the coupled 
equation of motion (13): 


(1 + u) 

m 


M „ F 2S i“i 

« 

0 

4 1 

V 

y 


IM L 0 1 

2C 2“2 U . 

X 2 



F o* (t) 


(B-l) 


An analog solution to Equation (B-l) above has been given in Appendix A. 
The general solution takes a complex form, with real and imaginary 
parts. One approach to obtaining the analytical solution to Equa- 
tion (B-l) is by first deriving and using the modal coordinates of the 
undamped system in Equation (B-l) to diagonalize the mass and stiffness 
matrices. This, however, will not diagonalize the damping matrix, but 
the common approximation of leaving out the off-diagonal terns will 
result in retaining only the real part of the solution. This is an 
acceptable approximation for small damping. 

Let us first write the frequency equation for the damped system 
since properties of the damped frequency equation will be used. The 
form of the solution is 



(B-2) 


Substituting in the left-hand side of Equation (B-l) results in the 
damped frequency equation: 

(^) * 2[R 'l * !l + l,)5 2 ] (0 * f" 2 * l ‘ H5 l S 2 * (1 * “’l fe) 2 


♦ 28^ ♦ « 2 R] * 


R = 0 (B-3) 


B-2 





The undamped frequency equation is then obtained by letting 


£ * £ ■ 0 : 
h *2 


a 1 * - [r 2 + (i + u)] ft 2 r 2 « o 


(B-4) 


where 


o * — 

W 2 


"i “i 

R = — = - 

W 2 ui ll + M 
s V s£ 


The roots to Equation (B-3) are w = S # and w = 0 w in which 11. 

ni 1 nd d d 1 

and 81,6 defined by 


_ _ A£1 

o, * n_ * — 

10 2 


AJ2 

0 * 0 + — 
2 0 2 


(B-5) 


and ii Q w 2 an d A0u)p, respectively, are the center and the beat frequency: 


w 2 / 2 

“> 2 % = -f ^ + 1) + P 


WgAQ = u /(R - l) + u 


(B-6) 


Associated with the frequencies w and obtained above are 
two undamped mode shapes, respectively defined by 




and {Yg} 


Next, the classical coordinate transformation is applied to x^ and x^ 
to diagonalize the matrices of Equation (B— 1 ) : 


*2} ■ m iy 


This results in: 



iL R a - f 

^ L 1 < 


+ { 2" 


X x [r (RC 2 ■ h 


A 2 |_r 


_1_ 

X 2 


(1 - fl 2 ) 2 

e i R T“ + ^ 


“nl ° Z 1 


0 <V Z 2 


where 


x = 1 [d - a ) + Mj 

q 

1 

x = [d - a 2 ) 2 + pJ 

a 

“2 


^ K ' 


p 0 (t) 


are the generalized mass. 


■-ssss* 
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It is .'een that the off-diagonai terns in the damping matrix will 
vanish when the damping ratio is proportional to the frequency ratio. 

No such limitation will be imposed, and therefore, t w^/uig. How- 

ever, the magnitude of the off-diagonal terms will be assumed small 
enough to be neglected in subsequent calculations. This is equivalent to 
neglecting the imaginary part of the solution, a commonly used approxi- 
mation that results in decoupling Equation (B-9), so that 


1 + 2£ .a) Z, + ta -Z. 
1 nl nl 1 nl 1 


Ll5 

a l\ 


i - n; 


V ‘> 


Z 2 + 2 «„2V Z 2 * \ 2 Z 2 - — 


snx„ 

2 2 


F M (t) 


(B-10) 


where now: 


2o C , 
nl nl 


(1 - V 2 .) 2 

V — ?- +e ^ 


2io £ 

2 “n2 5 n2 = X 


(1 - ^ )2 

X r ~ + 


Furthermore, from the theory of equation, the sum of the roots of the 
damped frequency equation. Equation (B-3), is equal to the coefficient 

o 

of ( p / u> 2 ) . Therefore: 


{ nl“nl * C „2"n2 = 2Bu 2 


(B-ll) 


0 = | U-jR + (1 + w) c 2 ] 


where 


The following approximate relationship will be useful: 


5 nl“nl ' 5 n2 w n2 


20u). 


(B-12) 


where 


t b e (< - i + - * fc 


R + (l + u) C, 


Equations (B-10) may now be solved for e prescribed force system 
F Q£ (t). The specific case when F^t) is an inpulse delta function of 
amplitude applied at t * 0 is discussed. By definition, the impulse 
Fo^(t) = response of an initial velocity v ''Tumerically 

F QJl = Vq^ for an impulse applied to a unit mass, as it is the case for 
Equation ( B— 9 ) - The homogeneous solution of Equation (B-10) is: 


-Cut -Cut 

_ _ . ^ nr nr . — . . _ ^ nr nr — . /_ , 

Z - A e sin u t + B e cos w t (B-13) 

r r nr r nr 


— / 2 

where u = u /I - C = damped natural frequency and A , B are 
nr nr nr r r 

constants (r * 1 for launch vehicle; r = 2 for the spacecraft). 

The initial conditions for the system at t = 0 are 


x^O) = x 2 (0) = 0 

x 1 (0) = v Ql and x 2 (0) * 0 (B-lU) 

By combining the transformation of (B-8) with the general 
solution of (B-13) and the initial conditions in (B-lU), we have 


B-6 
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^(t) 


(at 


v [ 4 {i - 

- q i> L ”»I 


” ^-i ) -5 ,lll t 

L. nim , ~ 
e sin m ,t 

nl 


n ?d - a?) - { 


2' t "n2 


ui „ 
n2 


W n2 t - 1 

Sln W n2 t I 


* 2 (t) 


R 2 v 


01 


K - nf) 




r L “m 


sm ui t 
nl 


(B-15) 


_1_ e - C n2“„ 2 * 

01 

n2 


sin <*> „t 
n2 


] 
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APPENDIX C 

RECOVERY OF LAUNCH VEHICLE MODAL PROPERTIES 






A. EQUATIONS OP MOTION 


In this appendix, equations are derived whereby the modes and modal 

responses of the launch vehicle with an interface loaded by only a 

rigid mass can be recovered from the given modes and modal responses of 

the same launch vehicle whose interface is loaded by an elastic and a 

rigid-mass simulation of the spacecraft. Let the subscripts S denote 

cantilevered spacecraft modes, and the subscripts C denote the composite 

launch vehicle free-free modes with an elastically loaded interface^ see 

Figure C-l. The motion of the total composite system of spacecraft and 

launch vehicle in Figure C-l subjected to an external force (F ) is 

J 

governed by the set of equations: 


{ V + l*e*cl { V + " <W { V (C_1) 


where q^ are available from the launch vehicle analysis. 

The modes of the launch vehicle with a rigid-mass can be computed 
from Equation (C-l) by subtracting the modal reactions of the elastic 
spacecraft at the interface. This is analogous to Equation (2) of the 
main text in which the sign of the spacecraft reaction force {F^} is 
negative since the spacecraft that was originally put on the launch 

vehicle is now being removed. In this case, the generalized coordinates 

« 

are denoted by {q } instead of {q }: 
c c 


Cq) + 

c 

From Equations 


K W c ]{ 4 c } + f“c 1 (q c } = [*cjl f V ' [ *ci H ‘V 

( 8 ) and ( H ) : 


(C-2) 


{F.} = [M ]{R ) ♦ [m ]{$•} 

1 11 1 J.o o 


{K,} = [*. ,]fq’} 


« M iil ■ [“hIres * L (n n' 


S = 1 


(C-3) 


j 


+ 
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Figure C-l. Composite Launch Vehicle With F,lastically 
Loaded Interface 


C-3 


Also, from Equation (3) of the main text, the spacecraft is governed by: 


(^1 ♦ [2« 5 » a H4jl ♦ - -(m si I(R t ) (c-M 


Equations (C-2, C-3, and C-M are now combined to give th- governing 
equations for the rigid-mass-loaded launch vehicle response: 


[ 


an - [♦ ci HM ii w ic i> 

-'♦cilKsl 

-[m .][♦. ] 

SI 1C J 

-hi 


111 } 


(Ji •{*“}«> 


c\j cj 
3 

0 

0 

2 

—CO 


s 


In a concise form. Equation (C-5) is written as: 


[M]{q’} + [51(4’} + [a. 2 ] {q'} = {F} 


and its undamped homogeneous part is written as 

[M] { q* > + {u> 2 } {q } = o 


(C-6) 


EIGENVAIAJE SOLUTION 


In general, if the system of composite modes, C, contained 
S-spacecraft elastic modes, and L-launch vehicle rigid-and-ela ;tic 
modes C = S + L, the undamped eigenvalue solution of Equation (C-6) will 
result in (L + ?S) eigenvalues « L , Q g2 and associated eigenvectors 

<V> <V> {V s 2 >- 

Of these, two S-roots will be repeated, * G g2 , and are asso- 
ciated with the elastic modes of the spacecraft. The remaining L-roots 
are associated with the rigid-mass-loaded launch vehicle free-free modes. 


C-2* 
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some of which will correspond to rigid body nodes with sero frequencies, 
and the remaining modes will correspond to elastic modes with non sero 
frequencies. The (L ♦ 2S) set of modes collectively referred to as [V] 

tv] - {V L j V a i i V i2 ] (C-7) 

will transform Equation (C-6) to the diagonal form 

+ [£2 2 HQ} « 0 (C-8) 

with 

{q* } * {V] {Q} (C-9) 

C. MODE SHAPES AMD ACCELERATIONS AT THE INTERFACE 

The acceleration {R^} at the interface can be expressed in either 
the {q* } or {Q} coordinates. Thus, 


1^) - [♦Jfq’) (C-10) 

or 

<y - i^hq) 

t 

where and [t|^] are the mode shapes at the interface in the (q } and 

the (Q) coordinates, respectively. 

Because of Equation (C-9), 

{R.} = r^Hq'} = [<#>.] [V]{^} 

= [h'.Hq} 

then 

[H' i ] = [4> i l [V] (C-ll) 
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D. MODAL ACCELERATION TIME HISTORY (Q(t)) 

The tine history of the modal acceleration corresponding to the 
newly obtained rigid-mass-loaded launch vehicle modes ( q) are obtained 
from Equation (C-9) so that 


(Q> - [V^lfq > ( C-12 ) 

The inverse fV in Equation (C-12) need not be performed. It can 
be found from the orthogonality condition that led to Equation (C-8): 

(v t HmHv1 . (U 


so that 


IV' 1 ) - [V T 1 [5) 


(c- 13) 


Therefore, 


(Q> = [V T ]fM]{q’} - ( C-lh ) 

However, {q} rather than {q* } are usually available from the launch 
vehicle analysis. Computing Cq* } as a transient solution of Equa- 
tion (C-5) is straightforward, but is both costly and unnecessary. For 
this reason, two approximate alternative solutions are shown next that 
give very acceptable accuracy at low cost . 


1. Approximate Solution 1 

Since the right-hard sides of Equations (C-l) and (C-5) are 
identical, they may be combined to yield the following exact expression: 


[M] 


•* I 





• • 

q c 


I“sK> + [2C s 0) s ]{q;} 

* . 


(s Vc I(i c - i;> * I“jH 


(C-l 5 ) 


c-6 







Tor the rigid body modes u • 0, the last tvo terms of Equation (C-15) 

C 

are identically zero. However* for the elastic modes* these last two 
terms are small because they express differences between two quantities 
of global charter having little dependence on the presence of the rela- 
tively small spacecraft. On account of this reasoning* the last two terms 
of Equation (C-15) can be neglected with little loss in accuracy. This 
implies : 


, • 
* q and q * q„ 


q * q 
c c 


Also* 


(C-l6) 


• ' . . . * 
q = q and q = q 
s s s 


Now, Equation (C-15) is approximate ly written as 


K 2 H<> + [ 2€ s “s H4 s } 


(C-17) 


which, when solved with Equation (C-lU), yields 


(Q) = [V 1 ] 




(C-l8) 


The modal coordinates (q } and {q } of the elastic simulation of the 

s s 

spacecraft may be obtained from an inexpensive transient solution since 
they are usually small (of,, the order of 3 to 6 modes). Tn this case, 
the lower partition cf Equation (C-5) can be solved for (q } and fq^} 
approximately by replacing {q c > by {q^}: 


(ql> * [2C u Hq’> + [«?](q*) = 


-[m .)[♦. J ( q (t)} 
si ic c 


( C— 19 ) 


C-7 


The approximation above is acceptable because fq > is affected only to 

the second order by (q c ) or (q c >, and also because there is no tuning 

between w and u> since u never contains ui exactly, 
sec s 

Thus, Equation (C-18) is alternatively expressed as: 


<5> * (v T i 



' ( V - 


( C-20 ) 


2. Approximate Solution 2 

Let the subscripts cR and cE, respectively, denote the subset of 
rigid body modes and the subset of elastic modes constituting the total 
number of the composite modes c. Thus c * cR + eE. This second 
approximation is based on the premise that, except for the rigid body 
modal accelerations {q cR } and {q^}, the acceleration time history of 
all the elastic composite modes are negligibly affected by the presence 
of the elastic simulation of the spacecraft. Therefore, for the elastic 
composite modes: 


* { %E (t)} 


( C-21 ) 


However, for the rigid body modes. Equation (C-l) gives 



} = [♦ 


cRj 


,( V 


( C-22 ) 


Also, Equation ( C— 2 ) for the rigid-mass-loaded launch vehicle gives 



} = [♦ 


cRj 


11 V 


tW“V 


(C-23) 
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Solving Equations (C-22) and (C-23) » and making use of (C-3), ve get 


(lil - 4 cR1 HM il Ht leB l) <„> * (5 oB ) ♦ [♦ cB1 ]t» ls ](5 3 ) 


* <«*> 

In the approximation of Equation (C-2l*)» {q } is confuted in the same 

s 

manner from Equation (C-I9) of the previous approximation. As such. 
Equations (C-2l) and ( C-2U ) along with Equation (C-19) provide an approxi 
mate alternative to computing all members of iq'}. Then again. 

Equation (C-l4) is used to compute {Q}. 
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APPENDIX D 

DISPLACEMENT SPECTRA FOR THE GENERALIZED 
RIGID BODY MODES 




In this appendix, an expression is derived for the generalized 
modal displacement spectra, Q gi , that should he used for the rigid body 
modes (P. ■ 0) in place of Equation (23) of the main text. 

From the lower partition of Equation (13), 


X„ + 2fU„X„ + «l“x “ - X 
2 72 2 2 2 1 


(D-l) 


where, according to Equation (14), 


s, - a * v i. 


(D-2) 


The homogeneous solution of Equation (D-l) having initial conditions 
x 2 (0) = x q and Xg(0) = x Q is: 


h,., 'W 

x 2 (t) = e 




+ ^2 Ui 2 X o'^ s * n w 2^ + X 0 COS u> 2^ | (D-3) 


The particular solution of Equation (D-l) is: 


x^(t) = 


i r - 

" /- x i 

* *' t=0 


(t) e 


-C-au(t-t) 


sin u 2 (t - t) dt (D-U) 


where 


m. - /l - f 
2 2 2 

Thus, the total solution of Equation (D-l) is 


x 2 (t) = x 2 + x^ 


(D-5) 


Now, the spacecraft modal response q (t) to rigid body launch vehicle 

s 

input (R = 0) is given by 


\ <t] TTmT *2 (t) 

s£ 


(D-6) 


D-2 




S8^ ^ 




Ot 


and the generalized modal displacement Q . for such cases is given by 

9 Jc 

the maximum of Equation (D-6). Therefore, 


I'm 


*si 


si 
1 + M 


si 


* 2 (t) 


max 


(D-T) 


Where again, R = 0 and x 0 (t) are defined by Equations (D-3) through 

(D-5). 

The homogeneous solution x^(t) can be ignored by starting to 
compute the forcing function well before the actual event, and only 
Equation ( D— It ) needs to be considered. Then Equation (D-7) becomes 


I'm 


„ si 

Q - — 

1 + M 


si. 


Xg(t) 


max 


Furthermore, by definition, |x^(t)| is the relative displacement 

max 

shook spectra of x (t) = (l + M ^) q^(t) 




+ (1 + V 

max Wg 


where is the acceleration shock spectra. Then 


Q si = 


/m 


si 2 
W 2 


l = 1 , 2 , •••, 6 


This is the same limit that Equation (23) takes for R -*■ 0 . Thus, 


J i Oi 


(ui« v n ») = 


= A P 


R-i-0 


(D-8) 


where . 
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